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Abstract 


Spectra  of  atmospheric  motions  as  a  function  of  observed 
frequency  may  depart  significantly  from  the  spectra  as  a 
function  of  intrinsic  frequency  due  to  a  nonzero  mean  wind.  In 
this  paperC-we^examine_*j  the  effects  of  Doppler  shifting  on  a 
model  spectrum  of  atmospheric  gravity  waves.  e*4«^>J:o  gain 
insight  into  the  effects  of  Doppler  shifting,  we  fvav€y derived 
analytic  solutions  by  approximation  of  the  intrinsic  frequency 
spectra  and  the  gravity  wave  dispersion  relation.  Our  results 
reveal  that  Doppler  shifting  can  have  major  effects  on  the 
observed  frequency  spectrum  of  both  horizontal  and  vertical 
gravity  wave  energy.  For  levels  of  Doppler  shifting 
representative  of  the  lower  and  middle  atmosphere,  possible 
effects  include  a  substantial  enhancement  of  horizontal  energy 
density  at  higher  observed  frequencies,  a  corresponding 
reduction  of  the  vertical  energy  density  at  higher 
frequencies,  and  a  significant  transfer  of  vertical  energy  to 
observed  frequencies  above  the  buoyancy  frequency.  The 
predicted  effects  are  found  to  be  consistent  with  some  of  the 
features  of  the  observed  frequency  spectra.  I  '  n  s ; 
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1 .  INTRODUCTION 

Internal  gravity  waves  are  known  to  contribute 
significantly  to  the  atmospheric  fluctuation  spectrum  at  a 
range  of  scales  and  frequencies.  In  addition,  they  are  now 
understood  to  have  a  number  of  important  effects  \n  the  lower 
and  middle  atmosphere.  Foremost  among  these,  and  the  subjects 
of  many  recent  studies,  are  the  induced  drag  and  turbulent 
diffusion  arising  from  the  saturation  of  upward  propagating 
gravity  waves  in  the  middle  atmosphere  (see  Fritts,  1984,  for 
a  review) .  More  recent  studies  suggest  that  saturation 
processes  may  also  be  important  in  the  lower  atmosphere. 
Palmer  et  al.  (1986)  have  suggested  that  gravity  wave  drag  may 
be  required  to  account  for  the  mean  structure  of  the 
tropospheric  jet.  At  the  same  time.  Smith  et  al .  (1987)  have 
provided  evidence  that  part  of  the  gravity  wave  spectrum  may 
be  saturated  throughout  the  atmosphere,  implying  a  possibility 
of  induced  drag  and  diffusion  at  all  heights. 

The  importance  of  gravity  waves,  both  in  accounting  for  a 
significant  portion  of  the  atmospheric  motion  spectrum  and  in 
driving  the  large-scale  circulation  and  thermal  structure  of 
the  lower  and  middle  atmosphere  has  stimulated  a  number  of 
studies  of  the  dominant  gravity  wave  scales,  characteristics, 
and  effects.  Some  of  these  have  focussed  on  the  determination 
of  individual  wave  parameters  (Smith  and  Fritts,  1983;  Meek  et 


al.,  1985a)  or  saturation  processes  (Fritts  and  Rastogi, 
1985),  while  others  have  addressed  gravity  wave  climatologies 
and  effects  in  a  statistical  manner  (Vincent  and  Reid,  1983; 
Smith  et  al.,  1985;  Meek  et  al.,  1985b;  Fukao  et  al.,  1985; 
Vincent  and  Fritts,  1987;  Fritts  and  Vincent,  1987).  While  the 
former  studies  are  needed  to  understand  the  detailed  processes 
by  which  gravity  waves  propagate,  dissipate,  and  contribute  to 
drag  and  diffusion,  statistical  studies  are  essential  in  order 
to  understand  the  average,  integrated  effects  of  gravity  waves 
in  the  lower  and  middle  atmosphere  and  to  parameterize  these 
effects  in  global  circulation  models.  Thus  it  is  important  to 
understand  any  effects  that  may  complicate  or  limit  the 
interpretation  of  gravity  wave  spectra. 

A  prime  example  of  such  a  complication  of  the  spectral 
approach  is  the  Doppler  shifting  of  gravity  wave  frequencies 
due  to  a  nonzero  mean  wind.  While  the  effects  of  Doppler 
shifting  should  not  cause  difficulties  in  case  studies  of 
gravity  wave  motions  and  effects,  they  pose  serious  problems 
for  spectral  studies.  Indeed,  neglect  of  the  effects  of 
Doppler  shifting  may  have  lead  to  erroneous  inferences 
concerning  the  nature  of  the  motion  spectrum  (Gage  and 
Nastrom,  1985 ) . 

VanZandt  (1982)  proposed  that  the  small-scale 
fluctuations  of  wind  in  the  lower  atmosphere  can  be  described 
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by  a  model  spectrum  of  gravity  wave  motions  very  similar  to 
the  Garrett  and.  Munk  (1972,  1975)  model  that  has  been  so 
successful  in  describing  the  corresponding  fluctuations  in  the 
ocean.  The  model  spectrum  is  formulated  in  terms  of  wavenumber 
and  frequency,  with  the  spectral  shapes  chosen  to  fit  the 
observations  as  well  as  possible,  while  remaining  consistent 
with  the  gravity  wave  dispersion  relation.  Thus,  the  model  was 
rather  empirical.  Recent  studies  of  the  saturation  of  gravity 
wave  spectra  have  removed  some  of  the  empiricism  by  providing 
theoretical  support  for  the  shape  of  the  vertical  wavenumber 
spectrum. 

There  remained,  however,  a  major  difficulty  in  applying 
the  existing  model  to  atmospheric  frequency  spectra  because  of 
its  dependence  on  intrinsic  frequency,  i.e.,  the  frequency 
measured  with  respect  to  the  local  mean  flow.  Unlike  the 
ocean,  where  the  observed  frequencies  are  very  nearly  equal  to 
the  intrinsic  frequencies  because  the  mean  current  is  usually 
much  smaller  than  the  intrinsic  phase  speeds  of  the  dominant 
gravity  waves,  the  mean  wind  in  the  atmosphere  is  often  larger 
than  the  intrinsic  phase  speeds.  When  this  happens,  Doppler - 
shifting  effects  can  be  substantial  and  the  intrinsic 
frequency  spectra  cannot  be  expected  to  apply. 

With  the  above  discussion  in  mind,  the  primary  goals  of 
the  present  study  are  1)  to  emphasize  the  role  of  Doppler 
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shifting  in  altering  the  observed  frequency  spectra  due  to 
gravity  wave  motions  and  2)  to  show  that  the  changes  arising 
as  a  result  of  Doppler  shifting  can  account  for  some  of  the 
observed  features  of  the  atmospheric  frequency  spectra.  To 
facilitate  our  discussion  of  the  physics  of  Doppler  shifting 
and  provide  a  discussion  that  is  easy  to  follow,  we  employ 
approximations  to  the  observed  motion  spectra  and  the  gravity 
wave  dispersion  relation  that  permit  us  to  obtain  analytic 
results.  We  believe  this  offers  significant  advantages  over 
the  treatment  of  Doppler  shifting  by  Scheffler  and  Liu  (1985, 
1986),  which  is  more  accurate  and  complete,  but  less  intuitive 
and  illustrative.  We  begin  in  Section  2  with  a  discrete 
intrinsic  frequency  spectrum  in  order  to  illustrate  the  manner 
in  which  both  discrete  and  continuous  wavenumber  spectra  are 
affected  by  Doppler  shifting.  The  resulting  frequency  spectrum 
is  seen  to  depend  fundamentally  on  the  distribution  of  energy 
with  horizontal  wavenumber  (or  phase  speed).  In  Section  3,  we 
introduce  a  continuous  intrinsic  frequency  spectrum.  The 
shapes  and  variability  with  anisotropy  and  degree  of  Doppler 
shifting  are  examined  in  detail  in  Section  4.  It  is  shown  that 
realistic  magnitudes  of  the  mean  wind  can  cause  the  shapes  of 
the  observed  spectra  to  be  quite  different  from  the  intrinsic 
spectra  and  that  the  spectral  changes  implied  by  Doppler 
shifting  are  consistent  with  some  of  the  observed  features  of 


the  frequency  spectra.  The  conclusions  of  this  study  are 
presented  in  Section  5. 

2.  DOPPLER  SHIFTING  OF  DISCRETE  INTRINSIC  FREQUENCY  SPECTRA 
The  purpose  of  this  section  is  to  demonstrate  the  effects 
of  Doppler  shifting  on  a  discrete  intrinsic  frequency  spectrum 
for  both  discrete  and  continuous  horizontal  wavenumber 
spectra.  This  will  provide  a  basis  for  understanding  the 
effects  of  Doppler  shifting  on  more  realistic  spectra,  which 
are  treated  in  the  following  sections. 

We  first  consider  the  normalized  energy  spectrum 
corresponding  to  a  monochromatic  wave  motion  at  intrinsic 
frequency  »'  -  «0'  and  horizontal  wavenumber  k  -  kQ.  This  is 
given  by 

E(k,u')  -  S<k-k0)S(w'-«o'>,  (1) 

where  the  S's  are  Dirac  delta  functions.  The  consequences  of 
Doppler  shifting  of  this  simple  spectrum  are,  of  course, 
trivial.  If  the  wave  motion  occurs  in  an  environment  with  a 

mean  wind  u  normal  to  the  horizontal  intrinsic  phase  velocity 

2 

of  the  wave,  of  magnitude  c0'  -k0w0'/k0  ,  the  frequency  seen 
by  an  observer  at  rest  remains  unaltered.  However,  if  the  mean 
motion  has  a  component  in  the  direction  of  wave  propagation, 
the  observed  ( Doppler - shi f ted )  frequency  is  given  by 

«  «  u>'  +  k*u  -  «'  +  k  u  cos^ 

o  o  ' 

7  j 

■ 


(2) 


-  Wq(1  +  u  cos#/c^) 

-  w£(l  ♦  0Qcos# ) , 

where  kQ  is  the  (positive)  magnitude  of  k,  $>  is  the  angle 

between  the  horizontal  components  of  It  and  u  and  |30  -  u/c0'  is 

« 

the  scaled  mean  wind.  The  fractional  frequency  shift  is  given 
by 

w-w' 

-  -  B  cos#  (3) 

wo 

The  quantity  0ocos#  thus  provides  a  convenient  measure  of  the 
effects  of  Doppler  shifting  on  the  frequency  spectra  of 
observed  wave  motions.  It  is  generalized  below  to  a  spectrum 
of  waves. 

In  the  case  of  equal  and  opposite  wave  motions  at 
intrinsic  frequency  «Q'  and  with  horizontal  wavenumbers  kQ  and 
-k0,  Doppler  shifting  results  in  two  observed  frequencies  up- 
and  down-shifted  from  the  intrinsic  frequency  by  equal 
amounts.  For  a  more  general  gravity  wave  spectrum  with  a  range 
of  horizontal  wavenumbers  and  directions  of  propagation, 
energy  can  be  Doppler  shifted  to  a  broad  range  of  observed 
frequencies  above  and  below  the  intrinsic  frequency,  as  shown 


below. 
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We  now  consider  a  superposition  of  waves  all  with  the 
same  intrinsic  frequency  w0'»  but  with  a  continuum  of 
horizontal  wavenumbers,  with  an  energy  density  given  by 

&(<*>'  )xa 

E(k,«')  -  - Zr—  ,  (4) 

(i  +  4 

where  x  «  k/k*  and  a,  t,  and  k*  are  parameters  that  must  be 
chosen  to  fit  the  observed  spectra,  subject  to  certain 
constraints.  Indeed,  Garrett  and  Munk  (1975)  and  VanZandt 
(1982)  fitted  the  observed  oceanic  and  atmospheric  spectra, 
respectively,  with  a  «  0  and  t  -  2.4  -  2.5.  At  sufficiently 
small  x,  the  logarithmic  slope  of  E  is  a,  and  at  sufficiently 
large  x,  it  is  -t.  The  slope  changes  from  a  to  -t  near  k  -  k*. 
Moreover,  with  reasonable  values  of  a  and  t,  most  of  the 
energy  in  the  spectrum  lies  within  a  factor  of  three  of  k*. 
Thus,  k*  plays  a  critical  role  in  the  description  of  the 
spectrum.  For  this  reason,  we  have  chosen  to  call  k*  the 
"characteristic  horizontal  wavenumber".  It  is,  of  course,  the 
same  as  the  "horizontal  wavenumber  bandwidth"  of  Garrett  and 
Munk  (1975). 

In  order. to  obtain  the  energy  distribution  as  a  function 
of  observed  frequency,  we  must  express  k/k*  in  terms  of 
intrinsic  and  observed  frequencies.  From  (2) 

k  -  ( «  )/u  cos* .  ( 5 ) 
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The  gravity  wave  dispersion  relation  for  hydrostatic,  two- 

2  2  2 

dimensional  motions  with  f  <<  w  <<  N  is  k  -  w'm/N,  where  f 
and  N  are  the  inertial  and  the  buoyancy  (Brunt-Vaisala) 
frequencies,  respectively,  and  m  is  the  vertical  wavenumber. 
Then  the  characteristic  horizontal  wavenumber  k*  can  be 
related  to  «'  and  the  characteristic  vertical  wavenumber  m* , 
which  Garrett  and  Munk  (1975)  and  VanZandt  (1982)  assume  is 
independent  of  «' ,  by 


u'm. 


N 


(6) 


Then 


*7 


W-M' 


N 


c; 


ucos#  w'm^ 


- 

UCOS+ 


»'  0COS^' 


(7) 


where  c*'  is  the  characteristic  intrinsic  horizontal  phase 
speed  of  the  wave  field, 


.  W'  N 

c  ■  ,  at  — 

^  *  k  m 


(8) 


and  where  &  is  the  appropriate  generalization  of  the  scaled 
mean  wind  to  a  spectrum  with  a  continuum  of  wavenumbers  and 
phase  speeds. 


Then  (4)  can  be  written  as  a  function  of  intrinsic  and 


observed  frequencies  as 


E(co,«'  ) 


E(k,«') 


6  ( w'  )■ 


[<*-<*'  )/cos$  )  (  gfa)* 
(@W'+(tt-U' )/cos$) 


(10) 


It  is  now  trivial,  because  of  the  delta  function,  to  integrate 
over  w' ,  with  the  result 

( («-«' )/cos*)a(3»' ) t_1 

E(«)  -  - 2 - hrr~-  <n> 

(8«£+(tt-fc>£)/COS*) 

This  solution  is  valid  for  all  Indeed,  if  the  azimuthal 

distribution  of  wave  energy  were  specified,  it  could  be 
integrated  numerically  over  ♦.  But  since  our  objective  in  this 
paper  is  to  gain  an  understanding  of  the  effects  of  Doppler 
shifting  by  means  of  approximate,  analytic  solutions,  we  will 
consider  only  the  extreme  effects  of  Doppler  shifting  by 
choosing  ♦  -  0°  and  180°  so  that  cos$  »  +  1 .  Then 


relation  and  a  saturation  spectrum  of  gravity  waves  at  large 
vertical  wavenumbers  (Smith  et  al . ,  1987),  and  consider  a  - 
•1,  0,  and  +1.  These  three  spectra,  normalized  to  the  same 
energy  density  at  large  k,  are  shown  in  Figure  1.  The 
resulting  observed  frequency  spectra  will  obvio.usly  depend 
both  on  the  spectral  slope  a  and  on  the  degree  of  Doppler 
shifting,  expressed  through  3. 

We  examine  this  dependence  by  differentiating  (12)  to 
reveal  the  locations  of  the  peaks  in  E+(«).  These  occur  at 


wo(1  -  a  >  0 


,  a  *  0. 


(13) 


Thus,  provided  that  a  >  0,  the  peaks  are  each  displaced  from 
the  intrinsic  frequency  «Q'  by  a  frequency  increment  3a«0'/t. 
If  a  *  0,  the  peak  is  unshifted,  with  only  a  spreading  of 
energy  about  w0'  occurring.  Note  that  these  results  assume 
that  all  values  of  w  from  -•  to  +•  are  allowed.  However, 
observed  spectra  do  not  distinguish  between  positive  and 
negative  frequencies.  Thus,  E(w<0)  must  be  folded  back  on 
E(«>0)  to  form  E ( | « | ) .  This  will  alter  the  observed  frequency 
spectrum  only  slightly  for  small  3a/t,  but  will  result  in 
significant  modifications  for  3a/t  >  1 .  To  avoid  ambiguity, 
the  frequency  spectra  obtained  here  and  in  the  following 
sections  will  be  presented  in  the  form  E(|w|). 


*  • 

•  • 

The  frequency  spectrum,  E(|<*>l)>  obtained  assuming  equal 
up-shifted  and  down-shifted  energies  for  the  spectral  slopes 
specified  above  and  for  0  -  1  and  3  are  shown  in  Figure  2  a 

and  b.  It  is  clear  from  these  spectra  that  large  energy 
densities  at  small  k  cause  the  energy  to  remain  near  w0' 
whereas  small  energies  allow  the  spectral  peaks  to  move  as  a 
result  of  Doppler  shifting.  As  will  be  seen  in  the  following 
section,  however,  the  energy  content  form  of  the  Doppler- 
shifted  frequency  spectrum,  mE(|«|),  shows  that  significant 
energy  can  be  Doppler  shifted  to  higher  frequencies  in  spite 
of  the  peak  in  E( |w| )  remaining  near  the  frequency  at  which 
E(«' )  is  maximum. 

We  have  provided  in  this  section  two  examples  of  the 
effects  of  Doppler  shifting  on  simple  intrinsic  frequency 
spectra  and  discussed  the  causes  of  shifts  of  the  peak  energy 
densities  when  these  occurred.  Our  intent  was  to  provide 
physical  insight  into  the  consequences  of  Doppler  shifting  of 
atmospheric  gravity  waves.  Hopefully,  this  will  assist  the 
reader  in  understanding  the  effects  of  Doppler  shifting  on 
more  complex  spectra  discussed  in  the  following  sections. 

3.  DOPPLER  SHIFTING  OF  CONTINUOUS  INTRINSIC  FREQUENCY  SPECTRA 

Our  intent  in  this  section  is  to  derive  approximate 
expressions  for  the  frequency  spectra  of  horizontal  and 


vertical  gravity  wave  motions  arising  through  Doppler  shifting 
of  the  intrinsic  frequency  spectra.  As  in  the  previous 
section,  we  will  assume  that  the  gravity  wave  motions  are 
hydrostatic  and  two-dimensional  in  order  to  obtain  analytic 
expressions  for  the  Doppler  -  shifted  spectra.  We  will  also 
again  assume  that  wave  motions  propagate  only  with  or  against 
the  mean  flow  so  that  cost  -  +  1. 

3.1  Up-  and  Down-Shifted  Energy  Densities 
We  begin  by  assuming  a  horizontal  spectral  energy  density 
for  intrinsic  frequencies  w'  between  f  and  N  given  by 

Eh( k )  -  E  ( p- 1 ) fp~ 1  (14) 

°  k*(l  +  k/k*)’1 

where  p  and  t  are  the  spectral  slopes.  This  is  the  same  form 
as  the  two-dimensional  energy  spectrum  used  by  Garrett  and 

Munk  (1975)  and  VanZandt  (1982)  except  for  the  neglect  of 

2  2  2 
factors  (1  +  (f/w)  ),  consistent  with  the  condition  <•>  >>  f 

leading  to  (6).  Here,  E0  -  E+  +  E.  is  the  total  energy 

(kinetic  plus  potential)  and  E+  and  E.  are  the  portions 

associated  with  gravity  waves  with  intrinsic  phase  speeds  c’  * 

u'/k  a  0  and  c'  <  0,  respectively. 

We  assume  that  this  spectrum  is  independent  of  the 

background  structure  of  wind  and  temperature.  Since  almost 

nothing  is  known  about  the  shape  of  the  intrinsic  frequency 


spectrum,  this  is  at  present  the  only  justifiable  assumption. 
It  is  also  simple  and  convenient.  But  it  is  unlikely  to  be 
true  in  general  because  the  propagation  of  gravity  waves  is 
strongly  affected  by  gradients,  particularly  vertical 
gradients,  of  wind  and  temperature.  It  would  be  a  good 
approximation  only  if  the  processes  that  tend  to  drive  the 
intrinsic  spectrum  toward  an  equilibrium  shape  were 
sufficiently  fast. 

The  corresponding  vertical  energy  density  is  obtained  by 
using  the  continuity  equation  and  the  dispersion  relation  for 
two-dimensional  gravity  waves  to  write 


w 


2 


(15) 


The  vertical  energy  density  is  then 


E  ( k , w' ) 


Eh(k,«'  )~ 
uz 


-  ^(p-l)fP-1  W'2'P 


Liill 


N 


k*  ( 1  +  k/kA)fc 


(16) 


These  spectra  can  be  expressed  in  terms  of  w  and  w'  by 
eliminating  k  using  Eq.  (2).  The  resulting  transformed  spectra 
are 


Ei,u(»,tt'  )  -  Ei,n<k'w'  > 


ak 

aw 


— s. 


.t-i-p 


(3w'+(w-w' )/cos$ ) 


(17) 


EI,II(W'W,)  •  “vfP(t-l)0 


t- 1  « 


(ew'+(w-w' )/cos* ) 


■A  18) 


where  «h  «  Nzav 


E+ (p-l)f*.  The  Doppler-shifted  spectra, 


E^(«a)  and  Ev(w),  can  then  be  evaluated  by  integrating  over  w' 

from  f  to  N.  As  before,  there  are  two  cases  that  must  be 

considered: 

cost  -  +1  and  w  *  w'  (Case  I)  (19a) 

and 

cos*  -  -1  and  w  <  w'  (Case  JI) .  (19b) 

In  both  cases  (3»'  +  (w  -  w' ) )  is  positive  so  that  (  )fc  is 
real.  The  different  ranges  of  integration  over  «'  for  each 
case  are  depicted  in  Figure  3.  For  u  *  f  or  u  *  N,  the  limits 
of  »'  are  f  and  N.  For  f  *  u  *  N,  however,  the  limits  are  f 

and  w  for  case  I  (w  *  u' )  and  u  and  N  for  case  II  («<«')•  As 

before,  since  the  observed  sp'ectra  are  in  terms  of  |o>|  ,  the 
spectrum  in  the  IIA-  range  must  be  folded  back  onto  the 
positive  axis  to  form  E(|«|).  Then  the  contributions  to  E( |u| ) 
in  each  range  of  |<«>|  are  given  by 


E  +  E 

I IA+  IIA- 

(20a) 

eib  +  eiib  +  eiia- 

(20b) 

IC  +  eiia-* 

(  20c) 

Because  both  the  integrands  and  the  limits  of  integration 
in  range  B  are  different  for  cases  I  and  II,  we  will  consider 
each  case  separately  from  this  point  on.  The  horizontal  and 


vertical  energy  densities  as  functions  of  observed  (Doppler- 
shifted)  frequency  w  for  case  I  can  be  written  as 

r«' 

Ej(«)  -  J^U  Ej (w, w' )d u' ,  (21) 

where  wu'  (u  for  "upper")  equals  N  in  range  C  and  w  in  range  B. 
Likewise,  for  case  II 

rN 

E_.(w)  -  £_.(«,»' )d«' ,  (22) 

J«'  A 

where  «i'  (1  for  "lower")  equals  <*>  in  range  B  and  f  in  range  A. 

In  order  to  evaluate  these  integrals  and  determine  the 
effects  of  Doppler  shifting  on  the  frequency  spectra,  we  must 
specify  the  spectral  slopes  p  and  t.  Recent  atmospheric 
observations  suggest  that  t  -  3,  due  to  a  saturation  spectrum 

of  gravity  waves  (see  Smith  et  al.,  1987).  Observed  frequency 
spectra  of  horizontal  velocities,  on  the  other  hand,  exhibit 
considerable  variability,  but  are  most  often  near  p  ■  5/3. 
Unfortunately,  this  choice  does  not  yield  integrals  that  can 
be  performed  analytically.  For  this  reason,  and  because  we  do 
not  expect  the  results  to  be  sensitive  to  small  variations  in 
p,  we  choose  p  -  2.  (It  may  be  noted  that  convenience  of 

integration  was  one  of  the  reasons  that  Garrett  and  Munk(1972) 
chose  p  -  2  (Nunk,  private  communication)!) 

With  the  selection  (p,t)  -  (2,3),  (21)  and  (22)  take  the 

form 

E(w)  -  «Jxm/( a+bx ) ^dx ,  (23) 
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with  m  -  0  and  2  for  the  horizontal  and  vertical  spectra, 
respectively.  The  solution  for  the  horizontal  energy  density 
(m  -  0)  is  quite  simple,  with  one  term  in  each  region  of  the 
form 


That  for 


I 

the  vertical 


dx 

( a+bx )  ^ 
energy 


1 

2b( a+bx ) 2  * 
density  (m  - 


(24) 

2)  has  three  terms 


and  is  given  by 


I, 


x2dx 


a+bx) 


1 

b^L 


2a 


!n( a+bx )  +  ^  - 


(25) 


2 (a+bx) 

In  particular,  it  is  now  clear  that  the  horizontal  and 
vertical  energies  will  not  be  simply  related  at  a  given 
observed  frequency  because  that  relation  is  based  on  the 
intrinsic  frequency  of  the  motion  which  is,  in  general, 
unknown. 

We  now  proceed  to  use  (17)  and  (18)  in  (21)  and  (22)  to 
obtain  the  expressions  for  up-  and  down-shifted  horizontal  and 
vertical  energy  densities.  Before  evaluating  the  integrals, 
however,  it  is  convenient  to  express  all  frequencies 

A 

nondimensionally  by  dividing  all  frequencies  by  f,  ie.,  x  - 
x/f . 


Each  region  with  w  £  «'  (IB  and  IC)  contributes  a  term  of 
the  form  (24)  to  the  horizontal  energy  density  and  a  term  of 
the  form  (25)  to  the  vertical  energy  density.  Expressed  in 
terms  of  nondimensional  frequencies,  these  energy  densities 
may  be  written 


18 


eJ(w)  -  «he2(W^-i) 


(  2a>+  (  0  - 1 )  ( «^+l )  ) 


(«+P-l)  (<»>+(  P -1  )«^)' 


(26) 


Ej<«)  -  2a  f 2 — — 
V  (P- 


q2  (»+(P* l)»n 

m  ^ - H 

-lrl  l  w+e-i  J 


AM  A  A  Am  A 

3«Z+4«( P - 1 )«^  3<»r  +  4c*>(  3  - 1 ) 


2<«+(P-l)w^) 


2  ( u>+P  - 1 ) ' 


(27) 


Similarly,  each  region  with  «<<»>'  (IIB,  IIA+,  and  IIA-) 


contributes  terms  of  the  form  (24)  and  (25)  to  the  horizontal 


and  vertical  energy  densities,  respectively.  These 


contributions,  expressed  in  terms  of  nondimensional 


frequencies,  may  be  written 


h  2  .  -  ( -2»+(p+l)(N+»! ) ) 

E  _(<•>)  ■  a.P  (N-«() — * - « s m - * — 

(-»+(P+l)w')  (-«+(P+l)N) 


(28) 


«?,(•)  -  2°vt2— ~ 

(#+!)■ 


l-»+(P+l  )w£, 


A  M  A  A  A  M  A  A 

3«  - 4«( P+1 )N  3u  -  4w(  P+1 )  wj 


2(*w+(P+l)N)  2( -w+(P+l)w' ) 


(29) 


We  are  now  in  a  position  to  determine  the  horizontal  and 


vertical  energy  densities  for  any  |w|  using  (20),  (26)  -  (29), 


A  A 

and  the  appropriate  value  of  <•>]/  or  wu'  for  each  contribution. 


3.2  The  Doppler-Shifted  Spectra 


y*i  jL&ry&J  J  av^.O  v  v  >1  v  j 


The  Doppler - shi f ted  horizontal  energy  spectrum 
obtained  by  substituting  Eqs.  (26)  and  (28) 

^  A 

appropriate  values  of  <•>].'  and  «u'  into  Eqs.  (20a, 

A  A 

with  -«  replaced  by  <•>  in  the  Ejja-  terms.  The 
expressions  are 

t 

A  f  ( -w+e+i )2< -w+(8+i jn)2 


+ 


E.  ,  - 

(N-l) 


(  2c>>+(  fl+1 )  (N+l ) ) 
(w+8+l)2(w+(3+l)N)2 
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Eg(f*|u|*N)  -  -^(M-l) (2m^(P‘1? ) 


2  2 

(«+e-i)zw  z 


(-«+(3+l)N)  « 


(W+P+1)>+(P+1)N)^ 


and 

Ep(  |  w  |  >N )  -  Ttg2(N-l)  ■«-(2{,>+iE\1HN+1?J 
c  (w+e-i)^(w+(3-i)N) 

x  <w-»-e+i)z(w+(e+i)N) 


h( l«l )  is 
with  the 
,  and  c), 
resulting 


(30) 


(31) 


(32) 


The  Doppler-shifted  vertical  energy  spectrum  Ev( |w| )  is 
obtained  in  the  same  way,  by  substituting  Eqs.  (27)  and  (29) 


into  Eqs.  (20a,  b,  and  c).  Because  the  expressions  in  Eqs. 
(27)  and  (29)  are  so  lengthy,  however,  we  express  these  energy 
densities  as 

A  A 

E^(|«|<f)  -  EjI(w;«£-l)  +  EjX  (  )  ,  (33) 

f 

Eg(f*|«|<N)  -  Ej(«;w^-w)  +  Ejj(w;w'-w)  +  Ejj ( ) , ( 34 ) 

and 

A  A  '  A 

e£(|«!>N)  -  Ej(w?«^-N)  +  Ejjf-ttjw'-l)  .  (35) 

As  a  check,  we  set  6  *  0  in  Eqs.  (30)  -  (35).  The 

resulting  spectra  should  be  identical  to  the  intrinsic 

frequency  spectra  and  are  given  by 

Vi  ( E  +E  )«  .  E  .  . 

E||(f*!w|*N)  -  --<*>  -  ,  (36) 

e£(N<|»| )  -  e£( |»|<f)  -  0  ,  (37) 

Eg(f<|«|*N)  -  (E+  +  EJf/N2  ,  (38) 

and 

E^(  |  u»|  <f  )  -  e£(N<|»|  )  -  0  .  (39) 

These  expressions  are  indeed  identical  with  the  intrinsic 
frequency  spectra. 

It  is  also  instructive  to  consider  E^(|w|)  when  6-1.  In 
the  middle  of  range  B,  where  f  <<  w  <<  N,  the  horizontal 
energy  density  reduces  to 
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Thus  the  up-shifted  portion  of  the  gravity  wave  spectrum  makes 
the  predominant  contribution  to  the  enhancement  of  Eb^( |«| ) 
relative  to  the  intrinsic  frequency  spectrum.  Likewise,  the 
horizontal  energy  densities  near  w  ~  0  and  w  >>  N  are 

h  E 

|«|~0)  (41) 

and 

h  ( E  +E  _ )  * 

Ep(|<*>|>>n)  *  (42) 

« 

Only  the  down- shifted  wave  energy  contributes  for  |«|  <  f 
while  the  up-  and  down-shifted  energies  contribute  nearly 
equally  for  |w|  >>  N.  Also,  the  energy  density  at  large  |w| 
decays  rapidly  with  increasing  |m| ,  as  expected.  These 
characteristics  of  the  Doppler  -  shifted  frequency  spectra,  as 
well  as  their  dependence  on  0  and  anisotropy,  are  addressed  in 
detail  in  the  following  section. 

4.  DISCUSSION  OF  DOPPLER- SHIFTED  SPECTRA 

We  present  in  this  section  a  variety  of  results  obtained 
from  the  horizontal  and  vertical  Doppler - shi f ted  spectra 
derived  in  the  preceding  section.  Our  objectives  here  are  to 
illustrate  the  various  effects  of  Doppler  shifting  on  the 


horizontal  and  vertical  frequency  spectra  as  a  function  of  the 
scaled  mean  wind  3  and  to  point  out  some  of  the  errors  in 
interpretation  that  may  occur  if  the  effects  of  Doppler 
shifting  are  not  properly  taken  into  account.  For  the  most 
part,  we  consider  only  components  of  the  gravity  wave  spectrum 
that  propagate  with  or  opposed  to  the  mean  flow,  recognizing 
that  the  components  propagating  perpendicular  to  u  are 
unaffected  by  Doppler  shifting.  Thus,  our  two-dimensional 
spectra  can  be  viewed  as  upper  limits  to  the  effects  of 
Doppler  shifting,  with  the  energy  density  for  orthogonal  wave 
motions  corresponding  to  the  case  3  *  0. 

In  order  to  compute  the  Doppler-shifted  energy  spectra  we 
must  specify  E+/E.,  N/f,  and  3.  In  most  cases  we  assume  that 
the  two-dimensional  gravity  wave  spectrum  is  symmetric  so  that 
E+/E.  -  1.  For  N/f  we  use  a  nominal  value  of  100,  which  is 
appropriate  in  the  troposphere  at  middle  and  high  latitudes. 
We  also  note  that  3  -  u/(N/m*)  ranges  from  0  to  a  maximum 
value  that  is  uncertain  principally  due  to  the  uncertainty  in 
the  value  of  m*.  With  the  smallest  values  that  have  been 
quoted,  0.5  -  1.0  xl0‘^  cyc/m  in  the  troposphere  and  lower 
stratosphere  (see  Smith  et  al.,  1987),  3  can  be  as  large  as 
**  10.  Therefore,  the  resulting  two-dimensional,  Doppler- 
shifted  horizontal  and  vertical  frequency  spectra  are  shown  in 
Figure  4  for  3*0,  0.8,  2,  5,  and  10. 


Several  features  of  the  Doppler  -  shifted  frequency  spectra 
are  immediately  apparent.  First,  the  peak  in  the  frequency 

A 

spectrum  of  horizontal  energy  (Figure  4a)  remains  near  w  *  1 
( <ji  -  f)  as  0  increases,  as  anticipated  in  Section  2  for  a 
horizontal  wavenumber  distribution  of  the  form,  used  here. 
Again,  this  is  a  result  of  the  finite  energy  density  at  very 
large  horizontal  scales  for  which  Doppler  shifting  has  a 
negligible  effect  on  the  observed  frequencies.  Also  evident  in 
Figure  4  are  significant  modifications  in  spectral  shape, 
including  changes  in  spectral  power  and  slope  at  all 
frequencies.  There  is  a  tendency  for  the  slope  of  the 
horizontal  frequency  spectrum  to  increase  (toward  positive 
values)  as  0  increases,  though  only  weakly  near  N. 
Accompanying  this  change  in  slope  is  an  increase  in  horizo?  tal 
energy  density  at  higher  frequencies  and  a  corresponding 
reduction  near  f  due  to  the  preferential  Doppler  shifting  of 
energy  to  higher  observed  frequencies.  In  contrast,  the 
frequency  distribution  of  vertical  motions  is  seen  to  undergo 
decreases  in  slope  (toward  negative  values)  and  in  power  at 
higher  frequencies.  This  occurs  because  the  majority  of  the 
vertical  energy  resides  at  high  intrinsic  frequencies  and  is 
easily  Doppler  shifted  to  observed  frequencies  with  |w|  >  N. 

To  illustrate  the  sensitivity  of  the  observed  energy 
density  to  small  0,  we  examine  the  form  of  Eb^( |«| )  for  0  and 

A  A 

1  <<  w  <<  N.  Also  assuming  E+  -  E.,  (30)  yields 


Eg(  |w|  )  * 


(43) 


-^(i+62/(i+e) ) . 

f  <l>  ' 

h  “? 

For  small  0,  Eb  (  | co |  )  increases  as  3  because  changes  due  to 
Doppler  shifting  are  largely  off-setting  for  the  up-  and  down¬ 
shifted  components.  For  larger  3,  the  increase  becomes  linear 
because  both  up-  and  down-shifted  components  now  transfer 
energy  primarily  to  higher  frequencies.  The  increases  for  this 
isotropic,  two-dimensional  gravity  wave  spectrum  are  0.5, 
1.33,  and  2.25  for  3  *  1,  2,  and  3.  The  fractional  increases 
would  be  only  half  as  large,  however,  had  we  assumed  equal 
gravity  wave  energy  propagating  zonally  and  meridionally. 

The  fractional  changes  in  the  horizontal  and  vertical 
energy  densities  are  shown  for  the  different  values  of  3  in 
Figure  5.  These  curves  imply  a  systematic  increase  in  the 


ratio  Eh( 

!«l )/ev< |tt| ) 

for 

|  <*>  j  *•  N  _  as 

3  increases.  This  ] 

reduction 

in  Ev( |»| ) 

and 

increase  in 

Eh( | w j  )  at  higher  1 

frequencies  constitutes  a  departure  from  the  ratio  of 
horizontal  to  vertical  energies  expected  for  gravity  waves  at 
high  intrinsic  frequencies.  Such  a  departure  could  lead  to 
incorrect  inferences  concerning  the  nature  of  the  motion 
spectrum  if  the  effects  of  Doppler  shifting  are  not  taken  into 
account.  It  should  be  noted  that  enhancements  in  the  ratio  of 
horizontal  to  vertical  energy  densities  at  higher  frequencies 
and  negative  slopes  of  the  frequency  spectrum  of  vertical 


f 

I 
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energy  have  been  observed  in  situations  of  nonzero  mean  winds 
(Gage  and  Nastrom,  1985;  Rottger , 1981 ) . 

A  more  insightful  presentation  of  spectral  results,  in 
many  respects,  is  the  energy  content  spectrum,  «E(w)  (see 
VanZandt,  1985).  The  Dopple r - shi f ted  horizontal  qnd  vertical 
frequency  spectra  for  8-0,  0.8,  2,  5,  and  10  are  presented 
in  this  form  in  Figure  6.  These  spectra  give  a  more  accurate 
indication  of  the  distribution  of  energy  density  by  accounting 
for  the  increase  in  frequency  intervals  with  w. 

Considering  first  the  horizontal  energy  content  spectra, 
«E^(  | ci> |  )  (Figure  6a),  we  see  that,  whereas  the  peak  in  E^(  |w|  ) 
remains  near  f,  the  frequency  at  which  the  maximum  energy 
occurs  shifts  from  f  at  0  «  0  to  ~  lOf  at  0  -  10.  This 
illustrates  much  more  clearly  the  shift  of  energy  toward 
higher  frequencies  due  to  Doppler  shifting.  It  is  also  clear 
that  a  minor,  though  not  insignificant,  fraction  of  the 
horizontal  energy  is  Doppler  -  shifted  to  |w|  >  N  for  large  0 
and  that  relatively  less  energy  appears  at  |«|  <  f  than 
indicated  by  the  horizontal  energy  spectra  in  Figure  5. 

T.ie  vertical  energy  content  spectra,  wEv(|w|),  shown  in 
Figure  6b  are  likewise  very  instructive.  First,  this  graph 
shows  that  very  little  energy  is  Doppler  shifted  to 
frequencies  less  than  f,  as  might  be  erroneously  inferred  from 
Figure  4b.  Instead,  it  is  clear  that  a  major  fraction  of  the 
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vertical  gravity  wave  energy  can  be  Doppler  shifted  to  w  >  N 
.at  larger  values  of  0.  While  the  peak  horizontal  energy  is 
Doppler  shifted  to  |»|  ~  lOf  at  0  •  10,  the  peak  vertical 
energy  density  is  shifted  to  |u|  ~  2N.  Thus  measurements  of 
vertical  gravity  wave  energy  with  a  Nyquist  frequency  near  N 
will  alias  or  exclude  a  significant  fraction  of  the  vertical 
energy  at  large  0.  As  discussed  earlier,  the  more  severe 
effects  of  Doppler  shifting  on  the  frequency  spectrum  of 
vertical  wave  energy  are  a  consequence  of  the  occurrence  of  a 
majority  of  the  vertical  energy  at  high  intrinsic  frequencies. 

We  now  examine  separately  the  up-shifted  and  down-shifted 
contributions  to  the  energy  content  spectra  shown  in  Figure  6. 
These  are  illustrated  for  0  «  2  and  5  and  compared  with  the 
intrinsic  (0  -  0)  spectra  in  Figures  7  and  8.  One  motivation 
for  considering  each  component  of  the  gravity  wave  spectrum 
separately  is  that  the  spectrum  is  often  highly  anisotropic 
(Vincent  and  Stubbs,  1977;  Hanson  et  al.,  1979;  Fritts  and 
Vincent,  1987;  Vincent  and  Fritts,  1987)  and  may,  on  occasion, 
be  best  characterized  by  only  the  up-,  down-,  or  un-shifted 
component. 

For  both  horizontal  and  vertical  energy  spectra,  the 
major  contribution  to  departures  from  the  intrinsic  frequency 
spectrum  is  the  up-shifted  portion  of  the  gravity  wave 
spectrum  because  the  fractional  change  in  |o>|  is  larger  for 


motions  with  «  >  w' .  For  example,  for  gravity  waves  with  the 


characteristic  horizontal  wavenumber  k*  -  we  see  from 

(2) ,  with  $  ■  2,  that 


while 
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Here  the  up-shifted  component  changes  by  a  factor  0+1  while 
the  down- shifted  component  changes  by  0-1.  The  frequency 
difference  is  enhanced  for  c  >  c*'  and  reduced  for  c  <  c*', 
with  the  difference  between  up-  and  down- shifted  components 
remaining  the  same  as  0  increases.  However,  the  relative 
frequency  change  for  up-shifted  and  down-shifted  components  is 
reduced  as  0  increases,  with  both  components  ultimately 
contributing  to  a  shift  toward  higher  Doppler -shifted 
frequencies. 

Finally,  we  present  in  Figure  9  the  area-preserving 
energy  content  spectra  for  0  <  0  <  100.  (values  greater  than  ~ 
10  are  unlikely  to  be  realized,  but  are  included  to  address 
the  limiting  behavior.)  These  spectra  show  quantitatively  the 
fraction  of  energy  for  a  symmetric  two-dimensional  gravity 
wave  spectrum  that  is  Doppler  shifted  to  frequencies  outside 
the  allowed  range  of  intrinsic  frequencies  ( |«|  <  f  and  |w|  > 
N).  For  0-2  and  5,  these  fractions  are  both  ~  0.20  for  the 


horizontal  energy  and  ~  0.30  and  0.50  for  the  vertical  energy. 
Once  again,  because  most  of  the  vertical  gravity  wave  energy 
lies  near  «'  «*  N,  estimates  of  the  vertical  gravity  wave 
energy  are  more  likely  to  be  biased  by  Doppler  -  shifting 
effects  than  are  estimates  of  the  horizontal  energy.  It  is 
also  evident  from  Figure  9  that  as  0  becomes  large  the  spectra 

for  w  >  N  evolve  toward  a  constant  shape,  which  can  be  shown 

2 

to  be  E0 ( «/0 )/( 1  +  w/0)  .  This  can  be  understood  by 
recognizing  that  for  0  >>  1,  u  >>  c*',  implying  that  the  wave 
motions  are  effectively  imbedded  in  a  high-speed  flow  with  the 
frequency  spectra  given  approximately  by  a  Taylor 
transformation  of  the  horizontal  wavenumber  spectrum.  This 
explanation  is  consistent  with  the  shift  of  the  peak  in  the 
horizontal  energy  density  to  |<o|  -  0f  for  0  >>  1.  A  similar 
shift  of  the  peak  vertical  energy  density  to  |w|  0N/5  is 
also  observed. 

5.  SUMMARY  AND  CONCLUSIONS 

We  have  presented  a  simple  discussion  of  the  consequences 
of  a  nonzero  mean  flow  for  the  observed  frequency  spectra  of 
gravity  waves  with  prescribed  intrinsic  frequency  and 
wavenumber  spectra.  In  order  to  illustrate  the  effects  of 
Doppler  shifting  as  simply  as  possible,  we  wanted  to  present 
analytic  solutions.  For  this  reason,  we  elected  to  consider 


idealized  gravity  wave  motions  and  spectra.  We  also  chose  to 
consider  only  gravity  waves  propagating  in  the  direction  of, 
or  opposed  to,  the  mean  flow  since  these  two  cases  establish 
the  limits  of  possible  Doppler-shifting  effects.  Finally,  we 
assumed  that  the  shape  of  the  intrinsic  spectrum  is 
independent  of  the  background  wind  and  temperature  structure. 
In  spite  of  these  approximations,  we  believe  that  the  results 
of  this  study  provide  a  number  of  insights  into  the  role  of  a 
nonzero  mean  wind  in  determining  the  observed  frequency 
spectra  of  horizontal  and  vertical  gravity  wave  energy  in  the 
atmosphere . 

To  examine  the  dependence  of  the  observed  frequency 
spectra  on  the  distribution  of  gravity  wave  energy  with 
horizontal  wavenumber  (or  phase  speed),  we  considered  in 
Section  2  a  distribution  with  intrinsic  frequency  of  the  form 
E(«')  *  For  wavenumber  spectra  with  energy 
localized  near  a  dominant  horizontal  wavenumber,  Doppler 
shifting  was  seen  to  result  in  frequency  spectra  with  two 
peaks  located  near  w  -  (1+S)w0'  resulting,  respectively,  from 
the  up-  and  down- shifted  components  of  the  gravity  wave 
spectrum.  For  gravity  wave  spectra  that  have  sufficient  energy 
density  at  small  k,  however,  both  up-  and  down-shifted  energy 
densities  are  spread  about  <*>  -  w0'  because  the  finite  energy 
at  very  small  k  remains  virtually  unaffected  by  Doppler 


shifting.  It  is  important  to  stress,  however,  that  the  energy 
content  forms  of  the  Doppler - shi f ted  frequency  spectra, 
wE( | to | ) ,  do  exhibit  an  increase  in  the  frequency  at  which  the 
maximum  energy  occurs  as  0  increases.  This  is  because  the 
energy  spectra,  E(|w|),  do  not  account  for  tha  increasing 
density  of  frequencies  at  larger  w. 

In  Section  3  we  considered  the  effects  of  Doppler 
shifting  on  an  energy  spectrum  similar  to  the  Garrett  and  Munk 
(1975)  spectrum,  which  has  been  shown  to  be  a  reasonable 
approximation  to  the  observed  atmospheric  spectrum  (VanZandt, 
1982).  The  results  were  seen  to  be  in  qualitative  agreement 
with  the  results  just  reviewed.  The  peaks  of  the  Doppler - 
shifted  frequency  spectra  of  horizontal  and  vertical  energy 
both  remained  near  the  peaks  of  the  intrinsic  frequency  energy 
densities  as  0  increased  due  to  the  assumed  constant  value  of 
E(k)  at  small  k.  But  like  the  results  with  more  simplified 
intrinsic  frequency  spectra,  the  peak  energy,  wE(|w|),  was 
seen  to  shift  toward  large  w  for  large  0,  with  the  peak  values 
occurring  near  w  -  0f  and  0N/5  for  the  horizontal  and  vertical 
energy  spectra,  respectively.  At  sufficiently  large  0,  the 
spectra  were  seen  to  be  given  approximately  by  a  Taylor 
transformation  of  the  horizontal  wavenumber  spectrum  by  the 
large  mean  wind. 


For  values  of  3  representative  of  the  lower  and  middle 
atmosphere  (3  5),  Doppler  shifting  was  found  to  cause 
significant  asymmetries  in  the  up-  and  down- shifted  components 
of  the  resulting  frequency  spectra.  Changes  induced  in  the  up- 
and  down- shifted  components  were  seen  to  be  largely  off¬ 
setting  for  3  £  1,  while  both  components  contributed  to 
changes  of  the  same  form  for  3  2. 

Because  the  shapes  of  the  horizontal  and  vertical 
intrinsic  frequency  spectra  are  quite  different,  the  effects 
of  Doppler  shifting  on  the  two  spectra  are  quite  different. 
For  |w|  *<*  N,  the  horizontal  energy  density  was  seen  to 
increase  due  to  a  transfer  of  energy  from  «'  ~  f  to  higher 
|w|  .  On  the  other  hand,  the  vertical  energy  density  was  found 
to  decrease  in  the  same  range  of  |w|  due  to  the  transfer  of 
energy  away  from  «*  «*  N  largely  to  |w|  >  N.  These  effects 
contributed  to  a  change  in  the  shape  of  both  horizontal  and 
vertical  energy  spectra  and  to  a  large  increase  in 
E^(  |«|  )/Ev(  |«|  )  near  |o>|  **  N,  consistent  with  frequency 
spectra  observed  in  the  atmosphere. 

Although  Doppler  shifting  must  have  effects  such  as  those 
described  in  this  paper,  the  shapes  of  E^( |«| )  and  Ev( |«| ) 
spectra  that  have  been  observed  when  u  (and  3)  are  appreciable 
are  usually  not  similar  to  the  model  spectra.  There  are 
several  reasons  for  this.  First,  the  observed  spectra  are 


contaminated  by  processes  other  than  gravity  waves, 

particularly  at  low  frequencies.  The  horizontal  spectra  are 

contaminated  by  planetary  wave  activity,  while  the  vertical 

spectra  may  be  contaminated  by  the  vertical  motions  in  slowly 

varying  mountain  lee  waves  or  potentially  by  the  projection  of 

« 

essentially  horizontal  motions  into  the  vertical  by  sloping 
isentropic  surfaces.  Second,  our  formulation  of  the  effects  of 
Doppler  shifting  is  highly  idealized.  Model  spectra  suitable 
for  comparison  with  observed  spectra  would  include  the  full 
gravity  wave  dispersion  relation,  allow  for  azimuthal  wave 
propagation,  and  permit  more  realistic  intrinsic  frequency 
spectra.  Indeed,  Scheffler  and  Liu  (1986)  have  developed  the 
formalism  for  such  calculations.  Such  a  study,  however,  is 
obviously  beyond  the  scope  and  intention  of  this  paper. 
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Figure  Captions 

Distributions  of  energy  with  horizontal  wavenumber,  E(k), 
for  low-wavenumber  slopes  of  a  =  -1,  0,  and  1. 
Distributions  of  energy  with  Doppler  -  shifted  frequency, 
E(|u|),  for  low-wavenumbe r  slopes  of  a  =  -1,  0,  and  1  and 
(a)  (3-1  and  (b)  3.  Note  the  separation  of  up-  and  down¬ 
shifted  energy  peaks  for  a  >  0. 

Regions  of  integration  in  the  (w,w'  )  plane  for  the  up- 
shifted  and  down-shifted  components  of  the  horizontal  and 
vertical  frequency  spectra. 

Horizontal  (aj  and  vertical  (b)  energy  densities,  E^( |u| ) 
and  Ev(|«|),  for  E+  -  E.  and  0-0,  0.8,  2,  5,  and  10. 
Energy  density  ratio,  Eh( | w| )/Ev( |w| ) ,  for  the  parameters 
in  Figure  4.  Note  the  large  enhancement  in  this  ratio  at 
higher  w  due  to  Doppler  shifting. 

Energy  content  spectra,  wE^( |u| )  and  «EV(|«1),  of  the 
horizontal  (a)  and  vertical  (b)  energy  densities  for  the 
parameters  in  Figure  4.  Note  the  shift  of  energy  to 
higher  frequencies. 

Horizontal  (a)  and  vertical  (b)  energy  content  spectra 
for  up-shifted,  down- shi  f. ted ,  and  symmetric  (2D)  gravity 
wave  spectra  for  0-2. 

As  in  Figure  7,  but  for  0-5. 


Area-preserving  energy  content  spectra  for  horizontal  (a) 
and  vertical  (b)  energies  for  the  parameters  in  Figure  4. 
At  large  0,  the  distributions  tend  toward  an  invariant 
quasi -Guassian  shape  that  can  be  interpreted  as  a  Taylor 
transformation  of  the  horizontal  wavenumber  spectrum. 
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